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ABSTRACT 

Intriguing sub-TeV tails in the pulsed 7-ray emission from the Crab pulsar 
have been recently discovered by the MAGIC and VERITAS Collaborations. 



^ \ They were not clearly predicted by any pulsar model. It is at present argued 

that this emission is produced by electrons in the Inverse Compton process 
. occurring either in the outer gap of the pulsar magnetosphere or in the pulsar 

*n ' wind region at some distance from the light cylinder. We analyse another 

00 \ scenario which is consistent with the basic features of this enigmatic emission. 



. It is proposed that this emission is caused by electrons accelerated very close 

■ to the light cylinder where the plasma can not saturate induced huge 



electric fields. Electrons reach energies sufficient for production of hard 7- 
ray spectra in the curvature radiation process. Due to different curvature 
, radii of the leading and trailing magnetic field lines, the 7-ray spectra from 

hj ■ 

■ separate pulses should extend to different maximum energies. The scenario 

can also explain the lower level 7-ray emission from the interpulse region 
(between PI and P2) observed in the Crab pulsar light curve. Moreover, we 
argue that pulsars with parameters close to the Vela pulsar should also show 
pulsed emission with the cut-off at clearly lower energies ('^50 GeV) than that 
observed in the case of the Crab pulsar. On the other hand, such tail emission 
is not expected in pulsars with parameters close to the Geminga pulsar. The 
model also predicts the tail 7-ray emission extending up to ^50 GeV from 
some millisecond pulsars with extreme parameters such as PSR J0218+4243 
and PSR J1823-3021A. 

Key words: stars: neutron — radiation mechanisms: non-thermal — gamma- 
rays: theory 
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2 W. Bednarek 

1 INTRODUCTION 

The modem satellite telescopes (Fermi-LAT, AGILE) have detected over a hundred pulsars 
at energies between ~100 MeV to ~ 10 GeV (e.g. Abdo et al. 2010a). These pulsars show a 
variety of light curves but their spectra are generally well described by a differential power law 
with an exponential cut-off at about a few GeV. In fact, such behaviour has been predicted 
by the models considered for the pulsed 7-ray emission, i.e. the polar cap (Ruderman & 
Sutherland 1975, Daugherty & Harding 1982, Venter et al. 2009), the outer gap (Cheng et 
al. 1986, Romani 1996) or the slot gap (Arons 1983, Harding & Muslinov 1998). However, 
observations of the Crab pulsar, performed with the MAGIC Cherenkov telescope, reported 
clear emission at energies ~25 GeV (Aliu et al. 2008) and the evidence of the signal at even 
larger energies (above ~ 60 GeV, Albert et al. 2008). This emission component, extending 
to a few tens of GeV, has been also confirmed by the Fermi-LAT observations of the Crab 
pulsar (Atwood et al. 2009). Farther observations with the MAGIC and VERITAS telescopes 
have measured the pulsed 7-ray emission from the Crab extending up to ~ 400 GeV (Aliu 
et al. 2011, Aleksic et al. 2011a), providing detailed spectral information concerning both 
7-ray pulses (Aleksic et al. 2011b). Whether this sub-TeV pulsed 7-ray emission connects 
naturally to the maximum at GeV energies or creates a separate spectral component is at 
present unclear. 

These results attracted much attention since they were not commonly expected in the 
standard versions of the pulsar models. More geometrically advanced versions of the outer 
gap model postulate the appearance of the power law tails in the 7-ray spectrum of the 
Crab pulsar caused by the comptonization of the synchrotron UV-X radiation from the 
gap by secondary cascade pairs (e.g. Hirotani 2011, Aleksic et al. 2011a, Lyutikov et 
al. 2011). However, detailed modelling of the pulse emission extending up to at least ~ 400 
GeV have not been provided up to now. It has been proposed that such sub-TeV pulsed 
7-ray emission should also appear as a result of comptonization of the soft magnetospheric 
radiation by pairs already in the pulsar wind beyond the light cylinder radius (Bogovalov 
& Aharonian 2000, Kirk et al. 2002, Aharonian et al. 2012). Lyutikov et al. (2011) have 
discussed possible production of sub-TeV 7-ray tails by electrons accelerated within the 
outer gap region in the curvature radiation mechanism, concluding that the process is not 
efficient enough for production of 7-rays above ~100 GeV if the outer gap electric potentials 
are considered. However, they also suggested that if the acceleration process is very efficient, 
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On the origin of sub-TeV gamma-ray pulsed emission from rotating neutron stars 3 

electrons might produce sub-TeV 7-rays in the curvature process. In this paper we investigate 
this suggestion in more detail considering a specific scenario for the origin of the hard 7-ray 
tail emission observed in the Crab pulsar. We propose that electrons can reach energies 
large enough for production of sub-TeV curvature 7-rays in the region very close to the 
light cylinder radius, where the induced electric fields can not be saturated by the pair 
plasma. We show that such scenario is consistent with the main features of the pulsed 7-ray 
emission from the Crab pulsar. Moreover, it gives clear predictions for possible appearance 
of such hard 7-ray emission in the case of other 7-ray pulsars. The basic features of this 
scenario can be tested by the future observations with the modern Cherenkov telescopes 
(HESS, MAGIC, VERITAS and also planned CTA). 

2 GENERAL SCENARIO FOR THE PULSED SUB-TEV 7-RAY EMISSION 

We assume that the pulsed structure of the 7-ray emission from rotating neutron stars 
(pulsars) is due to the emission from a single magnetic pole. For an illustration, we show in 
Fig. 1 two critical magnetic field lines (marked by LL - leading line and TL - trailing line), 
extending to the light cylinder and anchored in a single magnetic pole. These lines bound 
in the opposite direction. The dipole structure of the magnetosphere of rotating neutron 
star can extend only to the light cylinder radius, -Rlc = cP/2tt = 4.78 x 10''Pio cm, where 
P = lOPio nis is the period of rotation of the neutron star (NS) and c is the velocity of 
light. Close to (and beyond) the light cylinder, the structure of the magnetic field has to 
gradually change to toroidal one as described by the asymptotic solution of the oblique 
split monopole (Bogovalov 1999). The detailed structure of the transition region beyond the 
light cylinder is not well known. It is expected to be rather narrow with the characteristic 
thickness scaled with the light cylinder radius (see Sect. 2.1 in Petri 2011). Therefore, in 
our scenario we also scale the curvature radius of the magnetic field lines, Rc, at the light 
cylinder with the light cylinder radius. We introduce the scaling factor, rj, which links -Rc to 
i?LC according to Rc = tjRlc- Note that the leading and trailing magnetic field lines has to 
connect smoothly to the toroidal structure of the magnetic field in the wind region on this 
same distance scale in the transition region beyond the light cylinder (as e.g. modelled by 
Petri 2011). Therefore, the curvature of the LL lines have to be smaller than the curvature 
of the TL lines (see Fig. 1). Due to this difference the acceleration process of electrons along 
different lines is expected to be saturated at different maximum energies. We show that 
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4 W. Bednarek 

because of this difference in the curvature radii the energies of curvature photons produced 
by electrons accelerated in the region of the TL lines are expected to be significantly larger 
than those produced by electrons accelerated along the LL lines. Since the parameter rj can 
not be at present determined based on any theory, as an example we apply the values in 
the range 77 = 0.5 — 2, i.e. smaller value for the LL lines and larger value for the TL lines. 
The comparison of the 7-ray spectra calculated in terms of such scenario with the future 
observations of the cut-offs in the 7-ray spectrum of separate pulses of the Crab pulsar should 
constrain the value of the parameter rj, providing interesting constraints on the models of 
the pulsar magnetospheres close to the light cylinder. 

The considered scenario predicts that the peaks in the pulsar light curves should have 
different spectra and structure. We identify the emission from the LL lines with the pulse 
PI in the Crab pulsar light curve. On the other hand, the pulse produced along the TL lines 
is responsible for the emission of the pulse P2 in the Crab pulsar light curve. Acceleration 
of electrons can also occur along the intermediate magnetic field lines (IL lines) which are 
closer to the magnetic pole. In this case the acceleration process of electrons is expected not 
to be so efficient since the acceleration gap along IL lines is likely to be narrow. This region 
of the gap is expected to be responsible for the 7-ray emission observed from the interpulse 
region of the Crab pulsar light curve. 

In principle our scenario might be responsible not only for the hard tail 7-ray emission 
but also for the whole lower energy 7-ray emission observed from the Crab pulsar. In such 
a case, most of the 7-ray emission is expected to be synchronized in phase. On the other 
hand, a lower energy part of the spectrum (at and below an exponential cut-off) could be 
produced in another mechanism (operating closer to the pulsar) which had been already 
mentioned in the Introduction. In fact, some of the inner magnetosphere models (i.e. outer 
and slot gap models) predict 7-ray emission from the parts of the magnetosphere close to 
the critical magnetic field lines (LL and TL). These lines are directly linked to the emission 
region at the light cylinder considered in this scenario. Therefore, reasonable synchronization 
in phase, but not exact (see Abdo et al. 2010b, Aliu et al. 2011, Aleksic et al. 2011), between 
emission coming from these two different regions is expected. 
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Figure 1. Schematic representation of the pulsed 7-ray emission geometry in the case of a single magnetic pole model (not to 
scale). The magnetic field lines extend from the magnetic pole of the rotating neutron star (NS) up to the light cylinder radius 
(marked by the dotted curve) having a dipole structure. The direction of rotation of the NS is marked by an arrow. In such 
case, the leading magnetic field lines (marked by LL) have to change its geometry very significantly at the region of the light 
cylinder. On the other hand, the trailing magnetic field lines (TL) are stretched out. We also show the example intermediate 
line (marked by IL) which curvature is expected to be intermediate between those of the LL and TL lines, e.g its radius can be 
close to the light cylinder radius {rj = 1). We assume that electrons are accelerated efficiently close to the light cylinder radius 
within the so called light cylinder gap. The light cylinder gap is thick in the region of the LL and TL lines but it is thin in 
the region of the IL lines. The thickness of the gap is defined by the multiplicity of pair plasma arriving from the inner pulsar 
magnetosphere to the light cylinder radius. For LL and TL lines, the multiplicity is significantly below the critical multiplicity, 
fi << ficTit (see for details Sect. 3.1). For the IL lines, the multiplicity is fi > ficiit- Electrons, accelerated in the gap, suffer 
curvature energy losses producing 7-rays. 



We make use of the classical scenario for the pulsar magnetosphere defined in the work by 
Goldreich & Julian (1969, the Goldreich and Julian (GJ) model). In this model, the pulsar 
magnetosphere is filled with plasma which charge density is required to maintain co-rotation 
of the magnetic field lines. In fact, this classical scenario is expected not be viable in the 
part of the magnetosphere which is charge separated (in the region between the critical and 
the last closed magnetic field lines, see Goldreich & Julian 1969). However, it might be still 
valid in the rest of the open magnetosphere, i.e. along the magnetic field lines which are 
above the null surface (Spitkovsky 2011). Therefore, we limit our considerations to the part 
of the magnetosphere between the magnetic pole and the critical line. 

We assume that in the region close to the light cylinder a huge electric fields are generated 
in the magnetosphere rotating with the velocity close to the velocity of light. In a small 
region around -Rlc, these electric fields can not be saturated by the polarisation of the 
plasma since the required Goldreich and Julian plasma density should increase to infinity 
(see Goldreich & Julian 1969). The electric field intensity close to the light cylinder can be 
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6 W. Bednarek 
of the order of, 

Ei = ecB ^ 3005lc ~ 2.8 x IO^BuP^q^ eV cm-\ (1) 

where the magnetic field strength at the hght cyhnder is i?Lc = -Bns(-Rns/-Rlc)'^5 Rns = 10^ 
cm and -Bns = 10^^-Bi2 G are the radius and the surface magnetic field strength of the NS, 
e is the electron charge. 

On the other hand, electrons lose energy on the curvature radiation, 

= 2eV/(3i?c) ~ 4.2 X 10-^V/(^^io)^ eV cm-\ (2) 

where 7 is the Lorentz factor of electrons. By comparing energy gains of electrons with their 
energy losses on the curvature radiation {E{ = P^), we estimate the equilibrium Lorentz 
factor of electrons at the light cylinder gap, 

7eq = 9 X 10^(5i2r/VPio)^/^ (3) 

Electrons can reach such Lorentz factors provided that the width of the light cylinder gap 
is large enough. As we discuss below the gap along specific magnetic lines is constrained by 
the pair plasma production processes in the magnetosphere below the light cylinder radius. 

3.1 The thickness of the light cylinder gap 

The potential drop through the light cylinder gap has to be able to accelerate electrons to 
the equilibrium energies {E^q = m^'jeq)- As we mentioned above, the thickness of the gap 
depends on the density of pairs which enter the light cylinder region from the pulsar site. 
The calculations of the cascade processes in the inner magnetosphere of pulsars (below the 
light cylinder) show that the density of created pair plasma can overcome the Goldreich & 
Julian density by a factor of the order of /i = ri/(nGj,o = ^B/2txc) ~ lO^^^^-* (Harding & 
Muslinov 2011, Hibschman & Arons 2001, Ruderman & Sutherland 1975). These multiplic- 
ities are still about an order of magnitude lower than required in order to account for the 
observed radiation from e.g. the Crab Nebula (e.g. de Jager et al. 1996, Arons 2011). 

We assume that electrons are accelerated in the electric field of the order of that given by 
Eq. 1 on the distance scale corresponding to the thickness of the gap, Ai? = i?LC — -Rb, where 
i?b is the distance of the bottom border of the acceleration gap from the NS. The density 
of plasma arriving to the light cylinder from the inner magnetosphere, defined by /i, has to 
be lower than the Goldreich & Julian density at the light cylinder region, described by the 
full formula ncj = "n-GJ.o/il ~ {R/ RhcY] (see Goldreich & Julian 1969). Note that Goldreich 
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& Julian density diverges to infinity at the light cylinder radius. Thus, the acceleration 
gap has to appear even for very large multiplicity of pair plasma produced in the inner 
magnetosphere. Such scenario seems to be reasonable provided that the current induced by 
the dense plasma close to the light cylinder does not produce strong enough magnetic field 
which becomes comparable to the dipole magnetic field. Note that the GJ density formula 
has been derived based on the assumption of pure dipole structure of the magnetic field. 
It seems that this is still the case since in the GJ model the energy density in the plasma 
even at the light cylinder is smaller by a factor of ~10^ than energy density of the dipole 
magnetic field (see e.g. review by Spitkovsky 2011). 

To estimate the thickness of the light cylinder gap, we assume the validity of the formula 
for the Goldreich & Julian density and compare it with the density of pair plasma escaping 
from the inner magnetosphere (defined by coefficient /i), 

/X = [1 - (i?b/i?Lc)T\ (4) 

By inverting the above equation, we estimate the location of the bottom border of the gap 
and the thickness of the light cylinder gap, 

AR = R^c -Rb = RLcil- Vl - l/yu) = i?Lc/2/i. (5) 

We can estimate the critical value of the pair multiplicity for which the potential drop 
through the light cylinder gap is enough to accelerate electrons to the equilibrium ener- 
gies. The following condition has to be met, E{ ■ AR ^ E^q. The critical value of the pair 
multiphcity, derived from the above condition, is 

^ 1.5 X 10'P~J/'B%%-'/'. (6) 

For pair multiplicities lower than the above value, electrons can reach the equilibrium Lorentz 
factors (given by Eq. 1) and produce high energy 7-rays in the curvature mechanism. The 
pair plasma is not expected to be produced efficiently in the cascade processes in the inner 
pulsar magnetosphere along the last magnetic field lines, so called leading (LL) and trailing 
(TL) lines. Along these lines, the acceleration process of electrons and production of 7- 
rays are expected to occur already in the inner magnetosphere (e.g. the slot and outer gap 
models). In these regions the multiplicity of created pairs should be rather low in order to 
enable acceleration of electrons. On the other hand, the magnetic field lines closer to the axis 
of the magnetic pole (intermediate lines - IL) are expected to be more densely populated 
with the pair plasma (see e.g. the model by Ruderman & Sutherland 1975). Therefore, we 
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8 W. Bednarek 

conclude that the hght cyhnder gap in the region of LL and TL hnes is expected to be thick 
and accelerate electrons to the equilibrium Lorentz factors(7eq). On the other hand, the gap 
in the region of IL lines is expected to be significantly thinner. It can accelerate electrons 
marginally to 7eq or below this value. Electrons in this part of the light cylinder gap are 
not able to extract energy equally efficiently as those accelerated in the gap along the LL 
and TL lines. Thus, they produce 7-ray on clearly lower level as observed in the interpulse 
region of the Crab pulsar light curve. 

3.2 Gamma-ray production 

We assume that the magnetic field is significantly distorted from the dipole structure in 
order to met the toroidal structure in the pulsar wind zone. In such case, the curvature 
radius of the leading and trailing lines can differ from the light cylinder radius already in 
the light cylinder gap (see Fig. 1). We can estimate the characteristic energies of curvature 
7-rays produced by electrons which move with the equilibrium Lorentz factors through the 
gap (for the case of fi < /icrit), 

E; = 3/ic7eV(47ri?,) ^ 4507]'/^ Bli^/P'J^ GeV. (7) 

The curvature radius of the leading lines is expected to be smaller than i?Lc (i-e. V < ^) 
the trailing lines larger than i?LC (^ > !)• Therefore, we identify the 7-ray emission from the 
pulse PI in the Crab pulsar light curve as coming from the leading lines (lower maximum 
energies) and from the pulse P2 as coming from the trailing lines (larger maximum energies). 

The calculations of cascade processes in the inner pulsar magnetosphere envisage that the 
intermediate magnetic field lines are much densely populated by the pair plasma than the 
site lines along which gaps in the magnetosphere below the light cylinder radius can develop 
(e.g. Ruderman & Sutherland 1975). The pair multiplicities in the region of the leading and 
trailing lines are expected to be significantly lower in order to allow for the gap formation 
along these lines in the magnetosphere below the light cylinder (see the slot gap and outer 
gap models). We assume that along the leading and trailing lines the pair multiplicities are 
below yUcrit- Howevcr, along the intermediate lines (between the leading and trailing lines) 
the pair multiplicities are much larger, i.e. they are close to (or even above) Hcnt- Therefore, 
electric potential through the light cylinder gap along the intermediate lines may not be 
enough to accelerate leptons to the equilibrium Lorentz factors (7eq)- As a result, 7-rays are 
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expected to be produced less efficiently along the intermediate lines (the region between the 
leading and trailing lines) than along the leading and trailing lines. 

The observations of the giant radio pulses, and also detailed calculations of radiation 
processes in the inner magnetosphere (e.g. Timokhin 2010), argue for highly non- stationary 
production of pair plasma in the inner magnetosphere. Since the thickness of the gap along 
the intermediate lines is very sensitive on the multiplicity of pair plasma, it is expected that 
the acceleration of electrons in the light cylinder gap along the IL lines should produce vari- 
able 7-ray emission in the Crab pulsar interpulse region. Therefore, 7-ray emission from the 
IL lines should be anticorrelated with the appearance of the giant radio pulses. In contrust, 
the 7-ray emission from the pulse regions should not show any features of correlation with 
the giant radio pulses. 



3.3 Other constraints on the radiation scenario 

The characteristic distance scale (the mean free path) on which electrons lose energy on 
curvature radiation is, 

Ac = mec27eq/^'c(7eq) ~ 1-7 X lOy/^Pli/'/B^'' cm. (8) 

where 7eq is given by Eq. [3] and Pc is given by Eq. [2]applaying 7eq. Since Ac is much shorter 
than the curvature radius of the magnetic field lines, 7-rays from curvature process can be 
emitted only within a cone defined by the angle, 

A«cooi ~ Xc/vRlc ~ 3.6 X 10-''P^iy{7]'/^B%^) rad. (9) 

On the other hand, the angular spread of this emission is also limitted by the thickness of 
the gap, 

Aagap ^ Ai?/r7i?LC ^ l/(2/ir/) rad (10) 

For the parameters of known pulsars and the multiplicity of pair plasma arriving from the 
inner magnetosphere, of the order of 10^ — 10^, the angular spread of 7-ray emission is rather 
small, in agreement with the observations of the width of the pulses in the Crab pulsar light 
curve. We expect that the width of pulses should rather correspond to the angular extend 
of the magnetic field lines along which the 7-ray pulsed emission is produced. These lines 
should physically link the light cylinder gap to the gaps appearing in the magnetosphere 
below the light cylinder. Therefore, similar location of the 7-ray pulses in the pulsar light 
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curve is expected from the gaps below the hght cyhnder and those produced in the hght 
cyhnder gap. 

The coohng process of electrons due to curvature radiation becomes inefficient if their 
energies drop below certain value 7min- This value can be estimated by the comparison of the 
mean free path for electrons on curvature radiation with the characteristic distance scale for 
the origin of this radiation which is the curvature radius of the magnetic field, i.e. Ac ~ ^Rlc- 
This condition is satisfied for 

7min ^ 6.3 X lO^r/Pio, (H) 

Electrons with energies 7min escape from the light cylinder region into the pulsar wind. These 
electrons can be farther cooled by Inverse Compton scattering of soft radiation as considered 
by Bogovalov & Aharonian (2000), Kirk et al. (2002) and Aharonian et al. (2012). 

7-rays produced at the light cylinder might have large enough energies to be absorbed in 
the local magnetic field close to the light cylinder {'-fB — > e^). Such absorption process should 
be important when the mean free path for 7-ray photons (see Eq. 3.3a in Erber 1966) becomes 
comparable to the light cylinder radius. This condition is equivalent to the approximate 
relation, (i?^/2meC^)(-BLc/-Bcr) ~ 1/15. Note that the estimated number 1/15 (see also 
Ruderman & Sutherland 1975) does not differ significantly for the 7-ray pulsars in the range 
from the millisecond pulsars to the pulsars with the parameters (B,P) close to the Geminga 
pulsar. By inverting the above equation, we obtain the constraint on the surface magnetic 
field of the pulsars below which the curvature 7-rays can escape freely from the pulsar inner 
magnetosphere, 

B,, ^ I.OPT/v'^'. (12) 

This condition is not very restrictive for the classical and millisecond pulsars. 

By inverting Eq. 7, we obtain the constraints on the pulsar parameters for which curva- 
ture 7-rays can be produced with specific energies, 

B,, ^ 0.32P',/'EtS/v'^\ (13) 

where = lOOE'ioo GeV is the energy of the 7-ray photon. This formula indicates that hard 
tails in the pulsed 7-ray spectra at ~ 100 GeV can be only observed in the case of pulsars 
with rather extreme parameters such as observed in the Crab pulsar or some millisecond 
pulsars (for which P ~ 1 ms and B ~ 10^ G), provided that the curvature radius of the 
magnetic field at the light cylinder gap is not very far from the light cylinder radius. 
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4 PULSED 7-RAYS FROM THE CRAB 

Up to now, hard 7-ray tails in the pulsed spectrum of only the Crab pulsar has been dis- 
covered (Aliu et al. 2008, Aliu et al. 2011, Aleksic et al. 2011a,b). The emission from both 
pulses (PI and P2) differ significantly (Aleksic et al. 2011b). 7-ray emission above ~46 GeV 
is clearly weaker from the pulse PI in comparison to the pulse P2. Moreover, there are 
some evidences of the extended tail in the light curve after the pulse PI (Fig. 1 in Aleksic 
et. al. 2011b). The 7-ray spectrum of the pulse P2 extends up to ~400 GeV without any 
evidences of the cut-off. In contrast, 7-ray emission from the pulse PI is detected only to 
~200 GeV, with some evidences of the cut-off at higher energies. The width of the pulse PI 
is estimated on Aa^'"^ = 0.025 ± 0.008 and the pulse P2 on 0.053 ± 0.015 (Aleksic et al. 
2011a). 

We determine the basic parameters characterising our scenario by assuming the surface 
magnetic field of the Crab pulsar equal to -Bcrab = 4 x 10^^ G and its rotational period 
-Pcrab = 33 ms. For these parameters the equilibrium Lorentz factors of electrons accelerated 
in the light cylinder gap are 7^'^^'^ ~ 9.4 x lO^r^^^^. The critical multiplicity of pair plasma 
arriving from the inner magnetosphere is fJ^crif' ~ 525?7~^/^ and the characteristic energies of 
7-ray photons produced in the gap are _£;Crab ^ igQr/^/^ GeV. The scenario underlined above 
requires that the multiplicity of pair plasma along the leading (LL) and trailing (TL) lines 
is clearly lower than The widths of the pulses (PI and P2), Aa'^'^^^, should be lower 

than estimated in Eq. [TDl This last condition, Aa'^'^'^^ < Actgap, allows us to constrain the 
lower limit on the multiplicity of pairs, /i > l/2T]Aa''"'^^^, along the lines on which 7-rays are 
produced in the region of the pulses PI and P2. This lower limit is equal to fi > 20/ri for the 
pulse PI and fi > lO/rj for the pulse P2. For the above estimated parameters, electrons are 
accelerated within the gaps along LL and TL lines to the equilibrium Lorentz factors being 
in agreement with the angular extend of the 7-ray pulses from the Crab pulsar. Moreover, we 
conclude that the light cylinder gap along the intermediate (IL) lines has to be thin since the 
multiplicity of pairs along these lines is expected to be large (i.e. fi > yUCTit*^)- Therefore, the 
production of 7-rays along the IL lines in the region corresponding the interpulse region in the 
Crab pulsar light curve is expected to be less efficient in consistence with the observations. 

In order to better envisage the potential of our scenario, and convince that the expected 
7-ray emission features can be consistent with the observations of the Crab pulsar, we 
perform calculations of the curvature 7-ray spectra produced by electrons accelerated along 
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Figure 2. The example 7-ray spectral energy distribution from the leading lines (dashed curves), corresponding to the pulse 
PI in the Crab pulsar light curve, and from the trailing lines (solid curves) corresponding to the pulse P2. The spectra are 
calculated applying a simple model in which monoenergetic electrons with the equilibrium Lorentz factors (Eq. 3) produce 
curvature radiation in the magnetic field with the curvature radius defined by »; = 0.5 (for the leading lines) and 77 = 2 (for the 
trailing lines). Two pulsars are considered: the Crab pulsar (thick curves) and the Vela pulsar (thin curves). 

the LL and TL lines, corresponding to the region of the pulse PI and P2 in the Crab pulsar 
light curve. The pulsar electrodynamics close to the light cylinder is not precisely known. 
Therefore, the detailed calculations of the whole range of the 7-ray emission features are 
not possible on this stage. We consider a simple toy model which at least can envisage the 
possibility of production of curvature 7-rays with energies observed in the pulse regions of 
the Crab pulsar. We assume that electrons have monoenergetic energies (given by Eq. [3]) 
and the curvature of the magnetic field lines is described by the value of the parameter rj 
equal to 0.5 (for the LL lines) and equal to 2 (for the TL lines). The results are shown 
in Fig. 2. It is clear that for reasonable values of the curvature radii of the magnetic field 
lines, the curvature spectra can extend through the energy range observed by the MAGIC 
Collaboration from the Crab pulsar in the region of pulse PI and P2 (Aleksic et al. 2012). 

5 PREDICTIONS FOR OTHER PULSARS 

Let us consider in a more detail the prospects for detection of hard tails in the pulsed 
7-ray spectrum of other pulsars. According to our scenario, detection of the pulsed sub- 
TeV emission in the case of pulsars, with parameters less extreme than those of the Crab 
pulsar, is not so promissing. We consider two well known classical 7-ray pulsars, i.e. the Vela 
and Geminga pulsar (the strongest objects representing groups of pulsars detected in the 
GeV energies). The characteristic energies of 7-rays (Eq. [T]) estimated for the Vela pulsar 
(5ns = 3.4 X 10^2 G and P = 89 ms) and Geminga pulsar (1.6 x 10^^ p = 237 ms) 
are E^^^'^ ^ 2^7]^/"^ GeV and _E'Gemmga ^ Q.Sr^^/^ GeV, respectively. For the curvature radii 
suitable for the Crab pulsar (e.g. t] = 0.5 and 2), hard tail emission (above ~ 30 GeV) might 
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be observed in the case of the Vela pulsar (see Fig. 2) but not in the case of Geminga pulsar. 
In fact, There are some evidences of the existence of such hard emission in the pulsed 7-ray 
spectrum of the Vela pulsar already in the Fermi-LAT measurements (see Fig. 7 in Abdo et 
al. 2010c and Lyutikov 2012). Note also strong 7-ray emission from the interpulse region in 
the case of Vela pulsar light curve which suggest that pair multiplicity along the intermediate 
lines is below or not far from the critical multiplicity also along some intermediate lines when 
the production of pair plasma in the inner magnetosphere is not efficient (periods between 
giant radio pulses). The hard 7-ray tails extending above ~ 100 GeV might only appear 
in the case of pulsars with less extreme parameters than those observed in the Crab pulsar 
provided that the magnetic field lines at the light cylinder are characterised by more extreme 
values of the curvature radn. 

Interestingly, the model predicts that some millisecond pulsars might also show hard 
7-ray emission in the pulsed spectrum, similar to that observed in the Crab pulsar, since 
for them the characteristic energies of curvature photons are similar to that expected for 
the Vela pulsar, i.e. ^ MOr^i/^^g^VA^^^ GeV, where B = lO^Bg G and P = IPi ms. 
For example, we consider in a more detail two MSPs: PSR J0218-I-4243 and PSR J1823- 
3021A recently detected in the 7-rays by the Fermi-LAT telescope (Abdo et al. 2009, Abdo 
et al. 2011). PSR J0218+4243 (surface magnetic field 4x10^ G and the rotational period 
of 2.3 ms) is able to produce curvature 7-rays with energies f» ITr^^/^ GeV and PSR 
J1823-3021A (4.3 xlO^ G and 5.44 ms) E^ ^ 22r]^/^ GeV. Note that these values are similar 
to those derived for the Vela pulsar. With the magnetic field curvature at the light cylinder 
described by the parameters applied for the Crab pulsar, the 7-ray emission from the trailing 
lines of these two pulsars might be detectable by the present MAGIC stereo and the future 
Cherenkov telescope systems (HESS II, CTA) with the threshold at ~ 30 GeV. 



6 CONCLUSIONS 

We propose that electrons accelerated close to the light cylinder (in the so called light cylinder 
gap) can reach sufficient energies for the production of curvature 7-rays with energies up 
to a few hundred GeV in the case of pulsars with parameters similar to the Crab pulsar. 
The scenario is able to explain the spectral differences between the Pulse PI and P2 in the 
Crab pulsar light curve as a result of different curvature of the magnetic field lines at the 
light cylinder along the leading and trailing magnetic field lines. The broad gaps close to 
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the light cyhnder appear as a result of efficient acceleration of particles close to the light 
cylinder where the density of the pair plasma cannot saturate the electric fields. The model 
also predicts the appearance of lower level 7-ray emission in the interpulse region (betwen 
PI and P2) as due to the acceleration of electrons in the gap which is narrower along the 
intermediate lines than along the leading and trailing lines. The width of the gap along 
different lines is determined by the multiplicity of the pair plasma created within the inner 
pulsar magnetosphere along different magnetic field lines. 

One of the interesting features of the proposed scenario is the relation of the efficient 
7-ray emission to the magnetic field lines where the pair plasma density is low. Because 
of that, we predict that the hard 7-ray emission from the interpulse region cannot occur 
when the giant radio pulses are formed (caused by the non-stationary outflow of dense pair 
plasma through the pulsar magnetosphere). In fact, the giant radio pulses are expected to 
be anticorrelated with the level of the hard 7-ray emission from the interpulse region. 

The model also predicts hard pulsed 7-ray emission from pulsars with less extreme pa- 
rameters than that of the Crab pulsar. However, this emission is expected to extend to much 
lower energies. For example, in the case of pulsars with the parameters close to the Vela 
pulsar, and also some miUisecond pulsars (e.g. PSR J0218-F4243 and PSR J1823-3021A), 
the pulsed 7-ray emission is expected to extend only to ~ 50 GeV. These hard 7-ray "tails" 
should appear only in the region of the trailing magnetic field lines. Our scenario does not 
predict any hard 7-ray emission in the case of pulsars with the parameters close to the 
Geminga pulsar. Note that Inverse Compton pair cascade models for the sub-TeV 7-ray tail 
emission in the inner magnetosphere predict the cut-offs at much larger energies even in the 
case of the intermediate age pulsars since the primary cascade 7-rays are not expected to 
be absorbed so efficiently (e.g. Hirotani 2011, Du et al. 2012). This feature should allow to 
distinguish between these two models. 

In summary, our scenario predicts clear emission features for the pulsed 7-ray spectra 
from the classical pulsars with extreme parameters and from some millisecond pulsars. They 
are sometimes in contradiction to predictions of the popular Inverse Compton models for the 
origin of the sub-TeV pulsed 7-ray emission, in which 7-rays are produced by the secondary 
cascade pairs in terms of the outer gap scenario (Hirotani 2011, Lyutikov et al. 2011, Du 
et al. 2012) or by pairs outside the light cylinder radius in the pulsar wind (Bogovalov 
& Aharonian 2000, Kirk et al. 2002, Aharonian et al. 2012). The future observations of the 
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Crab pulsar (and also other pulsars) by the modern Cherenkov telescopes should allow us 
to distinguish between these different scenarios. 
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